PROCEEDINGS 

°F SCIENCE 



0^ 
O 
O 
(N 

a 



(N 



m 
(N 
cn 

o 

O 



A test of first order scaling in Nf =2 QCD: 
a progress report 



Guido Cossu* 

Scuola Normale Superiore, Dip. Fisica & INFN, Pisa, Italy 
E-mail: |g.cossu@sns.it 



Claudio Bonati 

Dip. Fisica & INFN, Pisa, Largo Pontecorvo 3, 1-56127 Pisa, Italy 



E-mail: bonati@df.unipi.it 



Massimo D'Elia 

Dip. Fisica & INFN, Genova, Via Dodecaneso 33, 1-16146 Genova, Italy 
E-mail: |delia@ge . inf n . it 



Adriano Di Giacomo 

Dip. Fisica & INFN, Pisa, Largo Pontecorvo 3, 1-56127 Pisa, Italy 
E-mail: [digiaco@df . unipi . it 



Claudio Pica 

Brookhaven National Laboratory, Physics Department, Upton, NY 11973-5000, USA 
E-mail: |pica@bnl . gov 



We present the status of our analysis on the order of the finite temperature transition in QCD with 
two flavors of degenerate fermions. Our new simulations on large lattices support the hypothesis 
of the first order nature of the transition, showing a preliminary two state signal. We will discuss 
the implications and the next steps in our analysis. 



The XXVI International Symposium on Lattice Field Theory 
July 14 - 19, 2008 
Williamsburg, Virginia, USA 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



A test of first order scaling in Nf =2 QCD: a progress report 



Guido Cossu 



1. Motivation 



It's a long standing issue whether the finite temperature transition in QCD (at zero chemical 
potential) from a confined to a deconfined phase is really a true transition or a simple crossover 
at finite fermion masses. This is not just an academic question if one wants to pursue the idea 
of a dual symmetry responsible for confinement, as suggested by experiments on the free quarks 
abundance in nature [jl]]. In this respect, the most natural hypothesis is the presence of a true phase 
transition even at finite masses rather than a crossover. 

For three or more degenerate fermions we know that a first order transition occurs in a region 
nearby the chiral point [§, ^]. In the 2+1 case there aie some hints in favor of a crossover scenario, 
at physical quark masses, from the study of the susceptibility of the chiral condensate ||]. The 
most debated case is QCD with two degenerate fermions: previous studies [^, 0,|8|] - performed on 
lattices with a modest extent - claimed that the chiral transition belongs to the 0(4) universality 
class, thus leading to a crossover at finite quark masses. Such a conclusion implies also the presence 
of a possible second order end point in the temperature-chemical potential plane that should have 
a clear experimental signature (see for example [^, Despite the efforts, no evidence for such 
an end point has been found so far at heavy ions colliders (BNL-RHIC, CERN-SPS, GSI-FAIR). 
For all the above reasons, to establish the order of the chiral transition in two-flavor QCD is a 
particularly important problem in lattice QCD which is still open and deserves a careful and deep 
analysis. Starting from our previous investigation, we report here about the progress made during 
the last year, using larger lattices. 



2. Current status 



This work is the last step of a long-term project. The first step [ ]11[ ] was a direct check of 
a second order scaling. Pisarski and Wilczek [ [l^ ] predicted, by means of a chiral model (thus 
assuming that only the chiral degrees of freedom are responsible for the order of the transition), 
that the QCD transition at the chiral point should be in the 0(4) universality class if a IR fixed point 
exists', i.e. if it is a second order transition, not excluding the first order case. 

A finite size scaling (FSS) analysis is the best way to address the problem of determining the 
critical exponents of a transition in a lattice simulation. The problem has two scales, the temper- 
ature and the bare quark mass, as shown e.g. in the equations of the scaling of the free energy 
density and its derivative, the specific heat: 

L/kT ~ Li<^{TL\''' ,amqLy") (2.1) 
C„ -Co ~ Lf^^^iTLl^^am^Ly") (2.2) 

Our strategy is to get rid of one of the two scales, namely the second one, by fixing its value 
amqL^/ = const, and looking at the FSS in the other variable. Obviously this strategy implies 
previous knowledge of so one could test consistency of data with a particular scaling hypothesis. 
The result of these consistency tests are shown in figure [T[ If the 0(4), or 0(2), scaling IS correct 
then all the curves should fall on top of each other. The conclusion based on these first observations 



'in the case of staggered fermions this universahty class reduces to 0(2) at finite lattice spacing. 
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Figure 1: Second order consistency check. Upper row refers to (9(4) scaling, lower to 0(2). Columns refer 
respectively to the specific heat and the chiral susceptibility from left to right. 



is that a second order universality class as predicted by chiral models is excluded. Notice also that 
the universality classes considered have a < 0, i.e. the specific heat should not grow with the 
volume, in contrast with data. Notice that 0(4) and 0(2) critical indexes are both very close to 
those of the U {2)i x U (2)r/U {2)v universality class, predicted as a possible alternative in case of 
a light rj' meson at the transition [|T3|]. Therefore our data exclude that universality class as well. 
The following step is to check directly the first order scaling for which some hints are found 



using approximate scaling laws in []1 Ip. The results of such a test ||14(], where j/, = 3, are shown in 
figure ^ The specific heat scales quite nicely with the first order hypothesis and also the chiral con- 
densate susceptibility is in agreement if one excludes the curve at = 16 arguing that it probably 
lies outside the scaling region since it has the smallest volume and the heaviest mass. 



3. First order scaling analysis 

The preceding observations were convincing enough to argue in favor of a first order transition. 
However several questions still remain: 

• where is the scaling with the volume of the peak at fixed mass, expected for some small 
masses? (here is difficult to quantify the term "small") 

• where are the double peaks expected in the observables histograms near the transition point? 
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Figure 2: First order consistency check of the specific heat (left) and the chiral susceptibility (right). 



Before looking at the simulations let's consider again the scaling laws at fixed bare quark mass and 
in particular equation If a second order transition is present at the chiral point then at finite 
mass everything is analytical and no divergence can arise, i.e. in the thermodynamical limit any 
dependence on L should vanish. We can expand in terms of the inverse of the second parameter 
and find that the leading term in the expansion of (^^ must be oc 1/ {amqU'')°'-l^^^'^\ In the case of 
a first order transition, where the equations are valid if the transition is really weak (as this is the 
case), a constant term (in volume) is present because of a peculiar cancellation occurring only in 
this case: 

Ci/-Co~flm^>i(TV) + V(^»2(TV) (3.1) 

the second term giving non-zero latent heat. 

The relative weight of the two terms is unknown a priori. It is perfectly possible that the 
singular, diverging, term is really small at the volumes explored. Nevertheless one should observe 
a shrinking with the volume in the width of the specific heat curve, in contrast with the crossover 
case where a constant curve is expected. 

In order to check volume effects we decided to dedicate a large part of our computational 
facilities to a run with standard staggered fermions at aniq = 0.01335 on a 48"^ x 4 lattice. This 
corresponds to a pion mass of about twice the physical value and to a spatial size of 13 — 14 fm. 
We present here our results at this stage. 

We simulated four different temperatures around the peak of the specific heat. In figure ^ we 
show the plaquette and the chiral condensate histograms. Data seem to present a weak signal of a 
double peak for betas 5.2719 and 5.2720. However statistics is still too low to draw any conclusion 
in this respect. A close look at the spatial plaquette susceptibility ('^ specific heat) (fig. ^ shows 
that the peak does not grow with the volume but the reweighted curve shrinks with the correct 
factor (see right figure). 

4. Conclusions 

The determination of the order of the Nf = 2 QCD transition is an interesting, rather tough 
problem, and we are still far from a conclusive answer. We reported our progress in understanding 
various aspects of the problem. First we concluded that a second order phase transition in the 
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Figure 3: Spatial plaquette (left) and chiral condensate (right) distribution histograms. 
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Figure 4: Spatial plaquette susceptibility (left). RescaUng of /3 axis (right) showing a good scaling of peak 
width. The — 16 lattice is out of the scaling region. 



0(4) or 0(2) universality classes at the chiral point has to be excluded. Secondly, investigating 
the possibility of a first order nature of the transition, we found several indications that could be 
the right answer. A direct consistency check of this hypothesis gave a positive answer, especially 
for the specific heat, an observable that does not imply any assumption on the symmetry behind 
the transition. Looking for metastabilities proved to be a more difficult task. Very weak signals 
of double peaks were found. If the transition is really first order this implies a very small latent 
heat, so that double peak structures are masked by simple thermal noise. A simple analysis of the 
scaling equations shows that in the first order case we should expect two contributions to scaling, 
one regular and one singular giving the latent heat. The actual data sets suggest that the regular 
term is much bigger than the singular one at the volumes and masses explored. We observe the 
expected shrinking of the specific heat curve with the exponents predicted for a first order. 

We need to increase statistics to have reliable ensemble with small errors and if our early 
conclusions are confirmed (eventually by simulations with an improved action and/or finer lattice 
spacing) the standard crossover scenario has to be changed. 
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This work was done using the apeNEXT faciUties in Rome during the time of several months. 
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